A unified treatment of current-induced instabilities on Si surfaces 
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We introduce a simple two-region model where the diffusion constant in a small region around 
each step on a vicinal surface can differ from that found on the terraces. Steady state results for this 
model provide a physically suggestive mapping onto kinetic coefficients in the conventional sharp- 
step model, with a negative coefficient arising from faster diffusion in the step region. A linear 
stability analysis of the resulting sharp-step model provides a unified and simple interpretation of 
many experimental results for current-induced step bunching and wandering instabilities on both 
Si(lll) and Si(OOl) surfaces. 



Si surfaces heated with a direct electric current un- 
dergo striking morphological changes induced by step 
bunching and step wandering instabilities. P, Q These 
phenomena are great interest not only because of possible 
applications for directed growth and nanofabrication, but 
also as physical examples of pattern formation on large 
(/Ltm) length scales in a system driven far from equilib- 
rium by a weak and externally controllable field. 

We present here a simple model that describes the 
interplay between the externally driven diffusion of 
adatoms by the electric field, and the intrinsic modu- 
lations and anisotropy of the diffusion pathways and the 
attachment kinetics on terraces and steps. It provides 
a new interpretation of the kinetic coefficients used in 
traditional step models and shows that negative kinetic 
coefficients can arise from faster diffusion in the region 
around a step. Many features of the instabilities seen on 
Si surfaces can be understood from this perspective. 

Step bunching is seen on vicinal Si(lll) surfaces when 
the current is properly directed normal to the steps. 
The uniform step train is initially stable when the cur- 
rent flows in the opposite direction. This instability has a 
mysterious dependence on temperature, with three tem- 
perature ranges between about 850°C and 1300°C where 
the stable and unstable directions are reversed. More- 
over, in temperature range II (about 1050°C to 1150°C) 
after heating for several hours with a current in the nom- 
inally stable (step-down) direction, the steps undergo 
a novel wandering instability with finite wavelength in- 
phase sinusoidal undulations in their positions. 0,l3[ 

Current-induced step bunching also occurs on Si(OOl) 
miscut along (110), but with some notable differences 
presumably arising from the alternating (1x2) and (2x1) 
dimer reconstructions on adjacent equilibrium terraces. 
At low temperatures step bunching is found for cur- 
rent normal to the steps in both directions, but involving 
paired double height steps. 0,0,0,01 

The behavior in the lowest temperature range I of 
Si(lll) is well described by a continuum diffusion model 
with nonequilibrium boundary conditions at sharp step 
edges if one assumes that adatoms acquire a small 



ositive effective charge e* and undergo biased diffusion 
from a field-dependent force F = e*E. However, it is 
not clear how to modify this picture to account for the 
reversals in the stable current direction at higher tem- 
peratures. Experiments 12] have ruled out the simplest 
explanation, a change of sign of the effective charge |^, 
and it seems likely that different boundary conditions are 
needed to describe the wandering instability. 

To gain some insight we introduce here a simple model 
that can give a more detailed description of processes oc- 
curring in the region around a step. The characteristic 
surface reconstructions seen on semiconductor surfaces 
strongly affect surface diffusion rates and pathways. We 
expect a different local reconstruction of bonds in the 
vicinity of a step. This suggests it could be profitable 
to view a step "dressed" by its local reconstruction as 
defining a region of finite width s (of a few atomic spac- 
ings a) where adatoms undergo effective diffusion with a 
diffusion constant Dg that can differ from Dt, the value 
found elsewhere on the terraces. |l3| For Si(lll) we can 
take Ds and Dt as isotropic and assume the step region 
has a fixed width s at a given temperature. 

Thus a uniform vicinal Si(lll) surface can be viewed 
as an array of two-region units, made up of the nth step 
region of width s and its neighboring lower terrace region, 
with width It ^ s. We assume that the straight steps 
extend along the y direction and the step index increases 
in the positive x (step down) direction, with x measured 
from the center of the step region. See Fig. (1). 

The biased diffusion flux of adatoms with density c 
takes the form: Jq, — —DaS7ca+Daica, where a = (t, s) 
indicates the terrace or step regions and f = Y /ksT. (We 
neglect evaporation and assume a constant positive effec- 
tive charge.) We first consider the steady-state solutions 
that arise when the electric field is directed normal to the 
steps and let / = f • i. We can ignore the small effects of 
step motion on the steady state adatom density field and 
determine c by simply requiring V- Jq = in each region, 
along with continuity of c and J at the fixed boundary 
at a; = s/2 between the step region and the lower terrace 
region. In almost all cases of physical interest, the field 
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FIG. 1: Geometry and density profile of the two-region model. 
The sloping solid lines in the upper part of the figure give 
the steady state density profile as a function of distance x 
from the center of the left step region, denoted by the verti- 
cal dashed lines. Shown is a highly exaggerated profile for a 
downhill force and slower diffusion in the step region, yielding 
step bunching in range I in Si(lll). Also illustrated with the 
dashed-dot line is the extrapolation of the terrace profile to 
the center of the step region, thus determining the parameter 
cf in Eq. Q. The lower part of the figure gives a side view 
of sharp equilibrium steps and their associated step regions. 



is sufficiently weak that fit and fs are much less than 
one, and the steady state profiles are piecewise linear. 
In particular the steady state terrace density is 
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for s/2 < X < I ~ s/2, with a similar linear expression for 
Cj(x). The (constant) adatom flux is given by 
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Here I = It + s \s the distance between the centers of 
adjacent step regions, R = Dt/Dg is a key dimension- 
less parameter that describes the relative diffusion rates 
of adatoms in the terrace and step regions, and Cg^ is 
the average equilibrium density when / = 0. This is also 
the density at the center of the step (and terrace) region, 
so one can view the step region as comprised of a clas- 
sical local equilibrium sharp step at x = surrounded 
by a symmetric local region of width s with a different 
diffusion constant, as illustrated in Fig. (1). 

Equation shows that the slope of the steady-state 
density profile is proportional to [R — 1)/. Thus, there 
are four general types of profiles, arising from a downhill 
force (/ > 0) or uphill force (/ < 0) combined with faster 
diffusion in the terrace region (i? > 1) or in the step re- 
gion [R < 1). We will use these results to make a precise 
connection between the physically suggestive parameters 
of the two-region model and appropriate boundary con- 



ditions in an equivalent symmetric sharp-step model giv- 
ing the same steady state terrace profiles. Both bunching 
and wandering instabilities can then be readily described 
using this standard framework. 

The general continuum boundary condition in the 
sharp-step model assumes small deviations from local 
equilibrium and introduces linear kinetic coefficients k± 
to relate c^ (or c^), the limiting lower (or upper) ter- 
race adatom density at the step edge, to the associated 
terrace adatom flux into the step: 



k±{ct 



(3) 



Here Ceq = c%[^ + Tk] with F the capillary length (pro- 
portional to the step stiffness) and k is the local step cur- 
vature. (31] For Si(lll) experiments are consistent with a 
symmetric model where ~ = k. 

In the usual quasi-stationary approximation, the dif- 
fusion field with boundary conditions given by Eq. Q is 
first calculated for fixed step positions, and then the step 
velocity is determined from the net local flux and mass 
conservation: 



w„Ac = h ■ [Jj - 3f] - drJs 



(4) 



Here Ac 



n- 



is the change in areal density 



when an atom joins the solid, Js = —DsdTCs + Ds{f-i)cs 
denotes a tangential or periphery diffusion flux along the 
step edge (the sharp-step analogue of parallel diffusion in 
the step region of the two- region model), and Cs — CeqS 
gives the effective number of ledge atoms per unit step 
length with diffusion controlled by Ds rather than Dt- 

It is natural to identify the terrace width in the appro- 
priate sharp-step model with I ^ It + s and to relate the 
limiting terrace density cf in Eq. (|3Jl to (0) , the extrap- 
olation of the two-region terrace profile in Eq. Q to the 
center of the step region, as shown in Fig. (1). Relating 
parameters in discrete and continuum models by extrap- 
olation is well known in other interface applications |l2l | . 

Using these results and the flux from Eq. Q in Eq. 
©, we find to lowest order in / our basic result: 



l)s. 



(5) 



This equation relates the fundamental parameters R and 
s of the simplest two-region model to the kinetic coeffi- 
cient k in an equivalent sharp-step model. |l3j 

Thus a positive kinetic coefficient k can arise from 
slower diffusion in the step region (i? > 1), in accord 
with the usual picture of an attachment barrier in range 
I. Indeed the extrapolated profile in Fig. (1) corresponds 
exactly to the linear steady state profile analyzed in 
if Eq. (O is used to relate parameters in the two models. 
As i? — > 1, we have k — > oo, and we arrive at the local 
equilibrium boundary condition with — Ceq. 

However, if diffusion is faster in the step region than in 
the terrace region [R < 1), we find a new regime with a 
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negative kinetic coefRcient. Equivalently, the character- 
istic length d = Dt/k is negative, but with d > —s/2. 

The possibihty of a negative kinetic coefRcient in the 
presence of a Schwoebel barrier was first suggested by 
Pohti and Villain, but with no derivation or dis- 
cussion of any physical consequences. We argue here 
that negative kinetic coefficients can play a key role 
in understanding current-induced instabilities on Si sur- 
faces. This has quite different consequences than a model 
with permeable steps. Consider a small perturbation 
5x„ (y, t) = Xn (y, t) - xl = ee'^*+*«2'+*"'^ + c.c. of the 
uniform step train. We report results for a linear sta- 
bility analysis of the sharp step model in the physically 
relevant limits of weak fields {fl <SC 1) and long wave- 
lengths [ql <C 1). An instability arises from a positive 
uj ^LOi (/, 0) -I- L02 (q, f, (f>), where 
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Step bunching is controlled by tJi . A pairing instability 
with maximum amplitude at cj) — t: is found for df > 0, 
or (R — l)f > 0, using the two-region model parameters. 
Thus, the profile illustrated in Fig. (1), produced by a 
step down current and slower diffusion in the step region, 
is unstable to step bunching, consistent with the usual 
interpretation .7] of range I in Si(lll). But a similar 
unstable profile arises from a step up current (/ < 0) 
along with a negative d oi k. 

UJ2 characterizes 2D step wandering. The first term 
in square brackets is always stabilizing and has its mini- 
mum value for in-phase wandering with (j) = Q. The next 
term, proportional to the field, has two contributions. 
The first, proportional to Dtc^gdfq^, describes a MuUins- 
Sekerka or Bales-Zangwill instability 0| induced by the 
terrace density field for df > 0. But as shown above, 
step bunching occurs under these same conditions. Wan- 
dering of the bunched steps is generally suppressed, as is 
seen experimentally in range I of Si(lll). 

However, the second contribution, proportional to 
DsCggSfq^, represents an alternate and quite general 
mechanism for step wandering that is operative when- 
ever there is a downhill force (/ > 0). Downhill step per- 
turbations are amplified by a field-driven downhill fiux of 
adatoms along the step edge with steady state density per 
unit step length approximated by = c^^s. Using Eq. 
(0) we see that the last term in Eq. is always positive 
even when d < and the MuUins-Sekerka contribution is 
stabilizing. 



Consider now the implications of these results for 
Si(lll). The low temperature experiments are well ex- 
plained by slower diffusion in the step region, consistent 
with the usual picture of an attachment barrier. 7| At 
higher temperatures in range II we suppose that changes 
in reconstruction could result in faster diffusion in the 
step region, implying a negative d or fc in the sharp- 
step model. Bunching then is predicted for / < 0, and 
step wandering for /> 0, in a gree ment with experiments 
and computer simulations |l7l| of such a model. The 
negative kinetic coefficient reverses the bunching direc- 
tion, which allows the general wandering instability from 
a downhill force to be easily seen. Indeed d < and 
/ > represents the only case where step wandering oc- 
curs with current in the opposite direction to that giving 
step bunching. Using parameters appropriate for Si(lll) 
we find that the most unstable wavelength is of order /im, 
in qualitative agreement with experiment. 

Since only relative diffusion rates are important, one 
could imagine the diffusion rates changing again at higher 
temperatures so that d > 0, possibly describing range 
III. In this scenario, the transitions between the differ- 
ent temperature ranges would be associated with local 
equilibrium behavior as R passes through unity, where 
no step bunching or wandering would be seen. 

These ideas also provide an interpretation of elec- 
tromigration results for the technologically important 
Si(OOl) surface miscut along (110). At equilibrium rather 
straight 5*^1 steps that run parallel to the dimer rows 
of the upper A terrace alternate with much rougher Sb 
steps that run perpendicular to the dimer rows of the 
upper B terrace [l6|. Moreover, diffusion parallel to the 
dimer rows is up to a thousand times faster at low tem- 
peratures 0. For driven diffusion with a current normal 
to the steps we thus have Df » Df, and we expect 
that this difference will dominate the physics of current- 
induced instabilities of Si(OOl). 

To apply our step-region ideas to this case, we imagine 
as before that a classical local equilibrium step resides in 
the center of each step region, but now let the downhill 
half-step region differ from the uphill half-step region, 
and assume that diffusion in each half-step region is sim- 
ilar to that in the nearest adjacent terrace. Defining i?* = 
D\/D\, with i = {A, B), we can let i?* differ from unity, 
thus generating asymmetric kinetic coefficients. We re- 
quire only that Df > Df, which seems quite reasonable 
since Df > Df. Thus Df - Df ^ ^f^{Df - Df) 
with ^f^ >0. Special cases of this model include clas- 
sical local equilibrium steps where R'^ = R^ = 1 and a 
symmetric step model where Df — Df. 

Experiments show that when a direct current is ap- 
plied to a configuration of alternating Sa and Sb steps, 
the steps move in opposite directions and step pairs form. 
With a downhill current one finds double height Db steps 
(consisting of an upper Sb step and a lower Sa step with 
a very narrow A terrace trapped in between) separated 
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by wide B terraces; the equivalent configuration with Da 
steps and wide A terraces is seen for an uphill current. 
On continued exposure to current, a step bunching in- 
stability of the double height steps is seen for current in 
either direction at low temperatures. 0, Q 

The initial step pairing can be most easily understood 
by calculating the velocity of steps in a configuration of 
equally spaced straight Sa and Sb steps. US] 

Using the 

flux given by Eq. ^ with the appropriate values of Dt 
and Eq. we find that the initial velocity of an Sb step 
can be written as vb = KAB{Df — Df)f, where Kab is 
positive and symmetric in A and B. The initial velocity 
of an Sa step is given by the same formula when A and B 
are swapped and hence has the opposite sign. Thus, for 
positive /, B terraces grow and A terraces shrink. Given 
the great difference in Df and Df, this process will con- 
tinue (as is shown by a general analysis with unequal 
terrace widths until the terraces have very different 
sizes, with the final width /' of the narrow A terrace in 
the Db step probably controlled by step repulsions 
(not taken into account in this version of our model). 

To explain the continued bunching of the double-height 
steps, we use the two-region model, with major terraces 
separated by a double-height step, which we treat as a 
single effective step region. However a minor terrace now 
resides in the center of the effective step region, and we 
must take this into account in our extrapolation analysis 
leading to Eq. jSJ. We can proceed as before if we note 
that the effective equilibrium density in the center 6°^ is 
linearly modified by the weak field from its value Cg^ at 
the "real" local equilibrium step near the lower boundary [10] 



of the effective step region, so that cj = 
s)/2]. Thus the analogue of Eq. ® is 13 



-eq 



[1 - f{l' + 



(8) 



where ' denotes an effective double step parameter. 

As shown above, with a step down current (/ > 0) 
B terraces grow and Db steps form. According to Eq. 
(0), step bunching occurs when df > 0. Thus continued 
bunching of the Db steps requires that in Eq. © is 
greater than zero, or > 2 + V / s. (Note that a local 
equilibrium assumption with R^ = 1 can give step 
pairing, but is inconsistent with further step bunching.) 
With step up current (/ < 0) A terraces grow and Da 
steps form. Continued bunching now requires that d"* is 
negative^ or R^ < 2 + I' /s. This inequality can be sat- 
isfied even if R^ > 1, and depending on the value of I', 
could hold under rather general conditions. At higher 
temperatures we expect values of R^ closer to unity due 
to thermal fluctuations. If R^ becomes less than 2 + V / s 
at some higher temperature, then Si(OOl) could exhibit 
behavior much like range II of Si(lll), with bunching 
only for a step up current, due to the (effective) nega- 



tive coefflcient for the Da step. Step wandering of the 
Db steps from a step down current, suppressed by the 
bunching at lower temperatures, would also be expected. 

Generalizations of these ideas and applications to ex- 
periments where the current is directed at an angle to 
the steps , along with comparison to results of Monte 
Carlo simulations will be presented elsewhere. 0| We 
are grateful to Ted Einstein, Oliver Pierre-Louis, and 
Ellen Williams for helpful discussions. This work was 
supported by the NSF MRSEC grant DMR 00-80008. 
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